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O ' Abstract 
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. We study a class of classical solutions of three-dimensional Af = 6 superconformal 

Chern-Simons theory coupled with U{N) x U{N) bi-fundamental matter fields. Especially, 
time evolutions of fuzzy spheres are discussed for both massless and massive cases. For 
the massive case, there are a variety of solutions having different behaviors according 
5^ I to the value of the mass. In addition to the time-dependent solutions, we analyze non- 

BPS static solutions which represent parallel M5/M5 or M5/anti-M5-branes suspended 
by multiple M2-branes. These solutions are similar to the fundamental strings connecting 
two parallel (anti) Dp-branes in perturbative string theory. A moving M5-brane and 
domain wall solutions with constant velocity that are obtained by the Lorentz boost of 
the known BPS solutions are briefly addressed. 
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1 Introduction 



Low energy world- volume theory of multiple M2-branes recently attracted much attention after 
the proposal of the Bagger-Lambert-Gustavsson (BLG) model [1,2]. In the BLG model, the 
theory consists of matter fields with Chern-Simons term of a novel 3-algebra gauge field. This 
model exhibits an explicit A/" = 8 supersymmetry with manifest 5*0(8) R-symmetry. Quanti- 
zation condition of the 3-algebra structure constant was pointed out in [1] which defines the 
level of the Chern-Simons term. The explicit Lagrangian was constructed by focusing on a 
3-algebra with structure constant f"-^'^'^ = e"'^'^'^ known as algebra. This theory was expected 
to describe two coincident M2-branes located on the "M-fold" [3] and various properties of this 
model have been explored [4]. One would expect that higher rank 3-algebras describe more 
number of M2-branes. However, under the assumption of positive definiteness of the metric, 
there is only one non-trivial example of the 3-algebra, just only A4, and its direct product are 
allowed [5,6]. In order to circumvent this no-go theorem, 3-algebras including Lorentzian met- 
ric [7-9], non-antisymmetric structure constant [10,11] were discussed. It was addressed that 
the ^4 algebra can be interpreted as 5*0 (4) gauge symmetry. This fact provides a valuable 
intuition for the ordinary gauge group description of the model. Indeed, in [12], the SO{A) BLG 
model were re-formulated as an SU{2) x SU (2) superconformal Chern-Simons- matter theory. 

Meanwhile, Aharony-Bergman-Jafferis-Maldacena recently proposed three dimensional A/" = 
6 superconformal Chern-Simons-matter theory as an alternative model of N coincident M2- 
branes in C^/Z^, orbifold (ABJM model) [13]. The gauge group is U{N) x U{N) and four com- 
plex scalar fields are introduced as the (anti) bi-fundamental representation of this gauge group. 
Two gauge fields and corresponding to each U{N)s have levels k and —k respectively 
and the Lagrangian of the ABJM model exhibits an SU (4) R-symmetry. For SU{2) x SU{2) 
gauge group, the A/" = 6 supersymmetry is enhanced to A/" = 8 and the ABJM model pre- 
cisely recovers the BLG model. It is also argued in [13] that the model is dual to M-theory on 
AdS^ X S'^/Zk at large- A^. A lot of works on this model have been studied [14]. 

On the other hand, as in the case of the world-volume description of D-brane configura- 
tions in perturbative string theory, the classical solutions in the M2-brane world- volume theory 
captures various properties of branes existing in M-theory. Remarkably, BPS solutions of the 
ABJM model have been analyzed by several authors [15, 16]. In these papers, a fuzzy funnel 
solution which preserves at least A/" = 3 supersymmetry was found. Similar to the Basu-Harvey 
analysis [17], the cross-section of this fuzzy funnel is fuzzy representing A^ M2-branes at- 
tached to an M5-brane. These BPS solutions represent static configurations of the branes and 
are unable to capture the dynamics of brane systems. In addition to these static solutions, it 
is interesting to study time-dependent and non-BPS configurations of branes in M-theory. 

Motivated by these facts, in this paper, we demonstrate that there are classical solutions 
in the ABJM model that represent time evolution of M2-branes or M2/M5-branes combined 
system. In order to study the time dependence of the M2-brane configuration, we examine 
its second differential equation of motion. In general, it is difficult to solve these non-linear 
second derivative differential equations, however, for some special situation, we can solve it 
analytically. In addition to these time-dependent solutions, we study static solutions that would 
be interpreted as M2-branes stretched between two M5s or M5/anti-M5-branes. Similar to the 
D3_LD1 system [18,19,21], we find that there are two types of solutions. One is a wormhole-like 
solution, the other is a cusp solution. Similar analysis were performed for D5_LD1 [22] and 
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D7±D1 [23] systems. 

The organization of this paper is as follows. In section 2, we briefly review the ABJM model 
and its half-BPS solutions. The equation of motion is also determined. In section 3, time- 
dependent oscillating fuzzy solutions are analyzed both for massless and massive cases. In 
section 4, non-BPS configurations representing multiple M2-branes suspended between parallel 
M5/anti-M5 and M5/M5 are discussed. Section 5 is our conclusion and discussion. Properties 
of the BPS matrices can be found in appendix A. In appendix B, moving fuzzy funnels and 
domain wall solutions with constant velocity are briefly discussed. In appendix C, a simple 
form of the time-dependent solution is presented. 

2 Equations of motion and BPS solutions 

In this section, we briefly review (2 + 1) dimensional A/" = 6 Chern-Simons theory proposed 
in [13]. This theory consists of U{N) x U{N) gauge fields with level k and —k coupled to 
(anti) bi-fundamental matter fields. We basically employ the notation of [20] but with different 
normalization of U{N) generators such that Tr(T"T'') = The bosonic part of the 

ABJM action is 

•S* = Skin + ^CS + Spot- (1) 

Here each parts are given by 

5kin = Jd^xTT[-iD^,Z^)iD'^Zy -iD^WA)iD^^WA)^, (2) 



Air 



An' y 
'S'pot = YY' I d X 



Af^d^Ax + '^A^A^Ax - A^d^Ax - jA^A^Ax 



(3) 



+ / d'x Tr 



{Z^Z\ + W^^Wa){Z''z\ - W^''Wb){Z^zI - W^^Wc) 

^{Z\Z^ + WaW^^^Z^Z"" - WbW^''){zIz^ - WcW^^) 
-2Z\{Z^Z% - W^^Wb)Z^{zI,Z^ - WcW^^) 

-2W^^{Z^bZ^ - WbW^^)Wa{Z^zI - W^^Wc) 

W^'^zIw^^WaZ^Wc - W^'^zIw^^WcZ^Wa 
-Z\W^^Zj.Z^WBZ^ - zW^^Z^Z^WrZ"^ 



(4) 



The world-volume metric is chosen such as rj^i, = diag(— 1, 1, 1). Here A^ and are U{N) x 
U{N) gauge fields, Z"^, W^'^{A = 1, 2) are bi-fundamental (N, N) representation of the U{N) x 
U{N) gauge group, k is an integer specifying the level of the Chern-Simons theory. The gauge 
covariant derivative is 

D^Z^ = d^Z^ + lA^Z^ - iZ^A^, 

D^Wa = O.Wa - iWaA^ + iA^Wa, (5) 
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and the field strength for is defined by 

Ff,, = d^A, - d,A^ + i[A^, AX (6) 

and similarly for A^. We normalized overall U{1) charges to be +1. This model exhibits a 
manifest SU{2) x SU{2) x U{1)r global symmetry which is in fact combined with the S'f/(2)^ 
and enhanced to S'f/(4)/j. For k > 2, this model has explicit A/" = 6 supersymmetry. The 
supersymmetry transformation of the component fields can be found in [15,24,25]. This model 
is expected to describe low energy effective theory of N coincident M2-branes in C^/Zfc. The 
world- volume coordinates {t, x^,x'^) are identified with the space-time coordinates {X^, X^, X"^). 

Once we drop the gauge fields and Wa, the BPS equation for = Z^{x2), which preserves 
a half of the original supersymmetry can be derived by performing the Bogomol'nyi completion 
[IgJeI or supersymmetry transformation [15] as 

a^Z^ = (Z^4Z^ - Z^4Z^) . (7) 

Let us consider an ansatz 

Z^ = f{x2)S^, (8) 
where are constant matrices satisfying the relation 

= S^Sl^S"^ - S^Sl^S^. (9) 

We call these S"^ "the BPS matrices". The matrices were first found in [26] and the explicit 
form are presented in the appendix A. The two complex scalar fields Z^ [A = 1,2) have 
the physical meaning of the transverse displacement of M2-branes along four-directions, say, 
(X'^, X"^, X^, X*) directions. If we assume real function / and using these BPS matrices, the 
BPS equation reduces to 

d2f{x2) = -y/'(^2). (10) 

A solution to this equation is easily found to be 




f{x2) = \ —ix2- Xo) 2,X2>Xo. (11) 



Because of the noncommutative property of Z^, the solution exhibits a fuzzy configuration. 
Indeed, this is just the fuzzy funnel solution with its cross-section fuzzy S^. This solution 
represents an M5-brane located at X2 = Xq extending along (X^, X^, X*^, X^, X*) and sharing 
the X^ direction with N M2-branes. 

On the other hand, to elaborate time-dependent or non-BPS solutions, we need to study 
second derivative equation of motion. The equation of motion for Z"^ with Wa = A^j, = A^ = 
is obtained as 

nZ^ = ^|3(Z^4)2Z^ + 3Z^(4Z^)2-2Z^4(Z^4)Z^ 

-2(Z^4)Z^(4Z^) - 2Z^(4Z^)4Z^} , (12) 



^In fact, the BPS equation presented here is coming from the D-term completion [16]. The F-term condition 
in our ansatz provides a trivial solution only. 
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where □ = d^d^. Note that this equation is Lorentz invariant. We can easily show that the 
BPS equation (I7j) is consistent with this equation of motion. Let us analyze a solution of the 
equation of motion with specific ansatz. As in the case of the BPS configurations, consider an 
ansatz 

Z^ = /(x)5^, f{x)eC (13) 

with the BPS matrices S^, then due to the property of the BPS matrices, the right hand side 
of the equation of motion becomes proportional to (see appendix) and the equation (|T^ 
reduces to the equation for /(x), 

Uf{x) = ^f{x)\f{x)\\ (14) 

On the other hand, it is possible to construct one parameter deformation of the ABJM 
model preserving an SU{2) x SU{2) x f/(l)i? global symmetry and M = Q supersymmetry. It 
is just the mass deformed ABJM model [24,26]. The equation of motion for the mass deformed 
ABJM model is 

□Z^ = ^|3(Z^4)2Z^ + 3Z^(4Z^)2-2Z^4(Z^4)Z^ 
-2(Z^4)Z^(4Z^) - 2Z^(4Z^)4Z^} 
Jjim I^B^^^A _ z^ZlZ^] + m^Z^, (15) 



where m is the mass of Z"^. Here, we have dropped the gauge fields and Wa- Plugging the 
ansatz f|T3|) back into the equation f|T5|) . we find the equation for the massive case, 

°/(^) = - + "^V(a^). (16) 

As in the massless case, the BPS matrices are linearized on both sides of the equation (1151) . On 
the other hand, the 1/2 BPS equation of the massive ABJM model is [16] 

d^Z^ = -'^ \^{Z''ZI)Z^ - Z^(4Z^)} + mZ^. (17) 

This is consistent with the massive equation of motion (1151) . If we assume that /(X2) is a real 
function, the BPS solutions of the equation ( fT6l) are given by 



/±(^2) = ^^(l±e— (18) 

where /+ is a domain wall interpolating a trivial vacuum and a non-trivial fuzzy vacuum 
while /_ is a deformed fuzzy funnel [16]. The ansatz (I13p preserves SU{2) symmetry. 

Note that any scalar fields satisfying the equations (fT2l) . (fT5l) should satisfy the "Gauss' 
law" constraints, 

(9^z^)4 - z^'id.zy = 0, 

(9^Z^)tZ^ - 4(9,,Z^) = 0. (19) 
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These constraints are automatically satisfied for the ansatz ffT^ provided 

f{x) = e''^g{x), g{x)eR, (20) 

with constant phase G R. 

A comment is in order. Because the equation of motions for massless and massive cases are 
manifestly Lorentz invariant, we can construct Lorentz boosted solutions from known static 
BPS solutions. A few examples of these boosted solutions are presented in the appendix B. 



3 Oscillating fuzzy 

In this section, we consider time evolutions of fuzzy in the ABJM model by solving the 
equation of motion with the ansatz of purely time-dependent function f{t). Several works on 
time evolutions of fuzzy spheres have been studied in the context of D-brane effective theories. 
In particular, collapsing fuzzy S"^ [27-29] and 5*^, [30] solutions were analyzed by the Dirac- 
Born-Infeld (DBI) equation of motion in the fully non-linear level. We will show that a similar 
collapsing fuzzy sphere exists in the massless ABJM model. Moreover, in the massive ABJM 
model, there are variety of solutions depending on the value of the mass parameter. 



3.1 Massless case 



Consider an ansatz = f{t)S^ with real f{t). The equation ( fTOl) reduces to 

-d^f{t) = 3af{t), 
where a = This can be rewritten as 

f\t) = -afit) + co, 



(21) 



(22) 



where cq is an integration constant and the dot stands for the time derivative. Assuming an 
initial condition f(t = to) = 0, f(t = to) = /o > 0, we have cq = a/g . The solution is valid only 
in the region < |/| < /o- The physical radius of the fuzzy 5*^ is 



N 



:Tr 



NT,. 



-Tr 



ct qA 



2{N 



-nt). 



(23) 



Here we have used the relation Tr[5']^5'"^] = N{N — 1) and physical coordinate = ^/^Z^ 



T2- 

[16]. The constant T2 is the tension of an M2-brane. A numerical plot of the solution to the 
equation fl?I]) is illustrated in fig.[Il This is an oscillating fuzzy sphere and we call this type 
of fuzzy sphere "collapsing" fuzzy sphere since it collapses to zero size then expands again. A 
similar analysis of collapsing fuzzy S'^ from the viewpoint of the DBI equation of motion in 
D2-brane world-volume was carried out in [27] where the time-dependent solution is described 
by Jacobi elliptic functions exhibiting a large-small duality of the fuzzy sphere radius. 

Let us evaluate the decay time of the fuzzy sphere by analytically solving the equation (!^ . 
This can be done by separating / > and / < cases. For / > case, the equation (1^ is 
integrated to give 



df 



± 



dt. 



(24) 



to 



5 




Figure 1: A numerical plot of the solution for the equation (^11) with k = 1, f(t = 0) = /o 
1, /'(t = 0) = 0. 



The left hand side of the above equation is evaluated as 



6' 2 



V/o''6'2 

Here we have used the definition of (incomplete) beta function 



/7/o < 1- 



(25) 



p1 poo poo 

B{a,h) = dt - tf'^ = rftt"-^(t + i)-"-^= / rft r"-^(t - 1)"-\ 

Jo Jo Ji 



B{z,a,b) = dt ^-'{1 -t) 
Jo 

< Re{z) < 1, a,b> 0. 



6-1 



dtf'-^t+l) 



-a—b 



1 

l-z 



dt f-^t-l] 



-a—b 



Thus the equation (HM is 

Here 1(2;, a, b) is the regularized beta function defined by 

B{z,a,b) 
I{z,a,b) = . 

B{a, b) 

From this expression, we find the following analytic solution in the region / > 



(26) 



(27) 



(2J 



f{t) = fo 



' 6' 2 



I2nfi t-tp 



(29) 



where / ^(2, a, b) is the inverse regularized beta function. The solution for / < region can be 
obtained by reversing the sign in front of /o in the equation (l29l) and changing t ^ —t. The 



^Thc full global solution can be obtained by taking the massless limit of the solution ([ST]) . 
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decay time T = td — to, f{td) = is evaluated as 



T 



6/0 V" 



B 



1 1 
''6' 2 



B 



1 1 
6' 2 



k r(i)r(i) 
12-/0^ r(|) 



^^^,,0^^ ~ 7.28595. 

r(i) 



(30) 



A conserved energy of this oscillating fuzzy 5*^ is also evaluated analytically as 

47r2 „A 47r2 



E = I d'xTr[S\S^] (^{djf 



fc2 



A;2 



■f^N{N-l) / rf^x 



(31) 



For = 1, the energy vanishes indicating the fact that this solution is purely due to the 
non-abelian nature of multiple M2-branes. 

3.2 Massive case 

Once we introduce a mass term for Z^, the ansatz = f(t)S'^ satisfies the equation 

127r2 



dtm 



-fit) 



P " k 

Assuming an initial condition /(to) = 0, /(to) = fo, the equation ( !32l) reduces to 

/' = -«(/'-/o)+/^(/'-/o)-7(/'-/o), 

where we have defined 



a 



1^ 



, /3 



Airm 



k 



-, 7 = m 



The equation of motion can be rewritten as 
-- -AaT{T - /( 



r 



a 



(32) 



(33) 



(34) 



(35) 



where we have defined t = p. We first consider / > region. In this case, the equation (|35|1 
becomes 



dr 



- fi) {r' + [fi - !) ^ + -Mf. - Pf^ + 7)} 



±2^/a [ dt. 

J to 



(36) 



The integral on the left hand side of the equation flHol) is separately evaluated for the following 
two cases, ( A) + (/^ - £) r + i {af^ - (3^ + 7) = has real solutions, (B) + (/^ - r + 
a ('^/o ~ /^/o + 7) = has no real solutions. For the case (A), /o should satisfy < /g < |^ 
while /o > |£ for the case (B). If /q = 0, there is only a trivial solution / = 0. 
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3.2.1 Case (A) 

In this situation, the integrand on the left hand side of the equation ( l36i) can be written as the 
following form 

(37) 



V-r{r-f^){T-g,){T-g,) 
Here, the constants gi,g2 are given by 



91 



92 



-fo+-]+fo\l-Spo 



a 



-/o+-)-/o\/-3/o^ + 



2P 
a 



a 



.2 , 2/5 



a 



(38) 
(39) 



The existence of the real roots gi,g2 is guaranteed by the constraint 

2/3 



< /o' < 



3a 



(40) 



In this region of /g, both gi and g2 are positive semi-definite. The behavior of the solution is 
different depending on the values of gi and g2- In this region, there are following possibilities 
according to the value of /q , 

(i) gi>g2> fl > 0, [o < fl < \^ 



(ii) 9i>92 = fl > 0, /o' 



(iii) 9i> fo > g2>0, 



6 a 



6 a 



<fl< 



2a 



(iv) 9i = fo > 92>0, fl 



1^ 
2a 



3 a 



(V) /o > ^1 > ^2 > 0, 

(vi) /o>^i = ^2>0, (^/o' 
Solutions for the cases (i)-(vi) are studied separately in below. 
• Case (i) gi > g2 > fo > 



For the case (i), the integration is evaluated as 
"^0 dr 



P V-i^-9i){r-g2){T- f^)T ^{gi - f^)g.. 



zF arcsin^ 



^92 fi-P 

fS92-P" 



Hi 



(41) 



fit) 

. 01/k 




Figure 2: An analytic plot of the solution 
chosen. 



m = 1, /o = 0.1, to = 0. Plus sign in ± is 



where 



91-92 fo 



9i - fo 92 ' 

Here F{ip, k) is an elliptic integral with modulus n defined by 



de 



--, F(arcsin(x), k) = sn ^x. 



F{lf,K) 

\ — sin^ 6 

Here sn(a;) is the Jacobi's elliptic function. From the expression fHTj) . we obtain 



sn^ 


=Fv^Vfi'2(5'i - 


-m-to)_ 


- 1 


/oSn2 


T\fa^/92{,9i 


-m-to) 


- 92 



(42) 



(43) 



(44) 



The solution is the positive root of the equation (1441) . For / < 0, the solution is given by the 
square root of (jH]) with minus sign and replacement t ^ — t. An analytic profile of the solution 
can be found in fig. El This solution represents a collapsing fuzzy sphere with the decay time 
T given by 



T 



(^V9i-fo92 ^2 



(45) 



Case(ii) g^>g, = f^>OJ^-l^ 



'0 — 6a 



For the case (ii), the equation of motion (l33l) is rewritten as 

1 



r 



108a 



{6ar-pyi3af'-2/3). 



(46) 



fit) 




Figure 3: An analytic plot of the solution (1491) with k 
g2 = 0.0233. The plus sign in ±-y/a is chosen. 



m 



1, /o = 0.3, to = 0,gi = 0.205, 



Because 3a p — 2(3 < 3gi — 2/? = at /q = the fact that the right hand side of the above 
equation should be positive semi-definite requires that the solution should be = /q = 
Although this is an extremum point of the potential, this point is a local maximum and the 
configuration is unstable. 



Case (iii) 



9i > fo> 92>0 



In the case (iii), the integration is evaluated as 



dr 



where 



arcsm ^ 



' 91-92 p-n 

91 - fo P -92'' 



ti2 



(47) 



9i>f> fl 



1^2 



1 9i - fo 92 



The solution is 



fit) 



9i - 92 fo ' 

92igi - /o)sn^ [±^/afo^gl - g2it - tp)] - figi - ^2) 
(^1 - /o)sn2 [±y^/o^^i - g2it - to)] - igi - ^2) 



(48) 



(49) 



The solution for / < region can be obtained as well as for the case (i). Note that when 
the numerator in ( H9l) goes to zero, the denominator also goes to zero and /^(t) remains finite 
value. This means that the solution does not collapse into zero size as can be seen in fig. [31 
oscillates between gi and /q meaning oscillating fuzzy 5*^ with finite radius. 



Case (iv) gi = f > ^2 > 0, f^ 



2 a 



In this case, the right hand side of the equation (133|) is positive semi-definite only at the 
point = /o = 1^. This is just the vacuum configuration found in [26]. 
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fit) 




Figure 4: An analytic plot of the solution fl52l) with k = m = 1, fo = 0.43, to = 0, gi = 0.128, 
(72 = 0.00517. The plus sign in ±y/a is chosen. 



* Case (v) > > > 0, i£ < < M 

In this case, the integration is evaluated as 
"^0 dr 



'p V-i^ - fo)i^ - 9i)ir - g2)r 
where 



V (/o - 92)91 



F arcsin^ 



/o > f> 9u «3 



fo-9i 92 
fo - 92 91 



The solution is 



fit) 



91 



ifo - ^i)sn2 ^v^^(/2 -^2)^i(i - ^0 



91 n - P 

P, - 91 P 



,f^3 ,(50) 



(51) 



(52) 



9i 



An analytic profile can be found in fig. HI Note that because the numerator in the right hand 
side of the equation (13^ is always positive, this solution does not collapse into zero size but 
oscillates within non-zero values of the radius. 



Case (vi) 



3 a 



In this case 91 = 92 



6 a 



I/q . The integration is 



dr 



fi \r - 9iW-'^i'^ - fo) 
Consider t > 91 case. Then the integration is evaluated as 

1 



^log 



(53) 



(54) 
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fit) 




Figure 5: An analytic profile for the solution (1331) with k = m = 1, Jq 



3tt 



tn = 0. 



where k = j^fo- From this expression, we have the following solution 



fit) 



fo 



;i + exp[±2v^(t-to)])^ 



4 1 — exp [±2y^cm(t — to)] + exp [±4y^cm(t — to 



r, it > to). 



This solution is neither oscillating nor collapsing but shrinking down to a finite radius at t 
Indeed, at large t, we have 

hm fit) = f . 

t— >oo 4 

An analytic profile can be found in fig.O 
3.2.2 Case (B) 

In this case, the integral on the left hand side in the equation ( l36i) can be rewritten as 

dr riB~f)/(A~f) A-B dt 



ifs JiB-mA-fi) ^/WM^it^-i,^)it^ + u^) 

where we have defined 



and 



i^ix) 
Fix) 

Gix) 



—X + f^x^ 



a 

A 
B 



X 



/3-2«/2 
1 

/3-2a/o2 



3if) - ^Kff) 



j(/o) + ^Hff) 
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The functions j{x), h{x) are given by 



j(x) = ax'^ — Px + j, 

h{x) = 3aV-5a/5x^ + 4a7x2 + 2/5V-3/57x + 72. 



(61) 
(62) 



F(B) „ G(B) 

^ ^ > 0, = > 0. 



FiA) 



G{A) 



(63) 



2a' 



Note that because /^^ — Aa'j = 0, j{x) is positive definite and due to the condition /q > 
we have P — 2a < 0. From these facts, we can see A > 0, B < and G{A) > 0, G{B) > 
0, F{B) < 0. Ahhough, there are two possibihties 



FiA) > 0, (0 < A < /o'), 
F(A)<0, (A>/2), 



(64) 
(65) 



we consider F{A) > case, otherwise the integral becomes complex valued. The integral fl37|l 
can be rewritten by the following change of variable. 



(66) 



Then the integral fl57p is 



A-B 



cdu 



A-B 



Here 



F(arcsin(A), K4) — F(arcsin(Ao), K4) . (67) 



2/iM 



A-P 



2\ 2 



/i^(l - A^) 



■\//i2 + Z/2|/i| |m| 



/i^(l-Ag). 



(68) 
(69) 



Note that ± in the expression of c is determined by whether t is positive (+) or negative (— ). 
Finding the explicit form of the solution /(t) from the equation ( 1671) is straightforward but we 
do not present it here. However, it is obvious that this solution is an oscillating fuzzy sphere 
which can be collapsed into zero size as we will see below. Summary of these solutions is found 
in table [H 

Because the nature of equations of motion comes from the 6th-order polynomial function 
for the scalar potential, the above characteristics of the solutions is inherently understood as 
the motion of a particle which obey classical equation of motion in the 6th-order polynomial 
potential. Setting off a particle with zero- velocity from a certain point that is given by /o, 
on the one dimensional potential curve, its trajectory shows the behavior of the corresponding 
solution (fig. EI). Since the potential blows up at infinity, it must be bounded. If the initial 
point is within the region (a), namely |/o| < \fpj2(x^ the particle begins to oscillate around 
the origin as in the case (i), except for the origin /q = 0. Then it passes across the origin in its 
cycle, the radius of the fuzzy S"^ collapses. In contrast, if the initial point is within the region 
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Table 1: Solutions for the massive case 




Figure 6: The behavior of the equation of motion is able to be understood as the behavior of 
particle obeying classical equation of motion with the 6th-order potential. 



(b) , the particle oscillates around one of the local minima at /o = ±A//9/2a, and thus the 
radius never collapses in this case, that is responsible for the case (iii) or (v). Within the region 

(c) , the amplitude of the oscillating motion is large enough to get over both local maxima at 
/o = ±>y /3/Qa, and the radius shrinks in the cycle again, responsible for the case (B). The 
shrinking radius solution of the case (vi) (/q = ±^3/5/2a) is rather particular. That is related 
to the climbing hill motion, starting from the point of the slope at the same level as the local 
maxima, to the local maximum of the same side, that needs infinite time to reach the summit. 
A simple form of the time-dependent solution which is valid for arbitrary values of /o is shown 
in AppendixO For both (A) and (B) cases, the energy is evaluated as 
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E 




6 



) 



2 



(70) 



One can easily show that this expression is positive definite for all the value of /o except the 
critical points /o = 0, ±J^aX where the solution is given by / = /o = const. 




In this section, we analyze purely spatial, static solutions of the ABJM model. These are 
generally non-BPS solutions but have interesting physical meaning. 

A similar analysis in D-brane systems gives a hint to our case. It is known that a BPS 
solution in the D-brane world-volume describes D-branes of lower dimensions. Especially, TV 
monopole solution in D3-brane world-volume can be understood as N D-strings attached to 
the D3-brane. In such D3_LD1 configuration, the dual description by D-string world-volume 
effective theory is also possible [21]. Multiple D-string world- volume effective theory is described 
by the non-abelian DBI action [31] and generically a solution exhibits fuzzy configurations due 
to its noncommutative matrix structure. The BPS solution of this D-strings effective theory is 
a fuzzy funnel whose cross-section is fuzzy S"^ and one divergent point along the D-string world- 
volume direction is interpreted as the location of a D3-brane that expands into the transverse 
direction of the D-strings. In addition to this BPS solution it is known that there are non-BPS 
double funnel solutions whose cross-section is fuzzy S*^. They have two divergent points within 
the finite segment along the D-string world-volume. These divergent points are interpreted as 
the location of two separated parallel 3-branes. According to how to choose the integration 
constant, there are two types of double funnel solutions. One is the wormhole-like solution 
describing D3 and anti-D3 connected by D-stings [18]. The other is the cusp solution studied 
in [19] representing A^ D-strings suspended by two D3-branes. For large- A^, the fuzziness of the 
configurations is effectively smoothed out and the energy and the R-R charge corresponding to 
the configurations precisely agrees both from the D3 and Dl point of views. 

Analogous configurations are expected to exist also in eleven dimensional M-theory. Indeed, 
in the following, we demonstrate that such solutions actually exist in the ABJM model. 

4.1 M5/anti-M5 solution 

Consider the ansatz = f{x2)S^ with a real function /. The massless equation of motion 




4 Non-BPS solutions 



reduces to 



dlf{x2) = 3af{x2) 



47r2 



(71) 



a = 



15 



At fixed x^, the configuration is interpreted as a fuzzy S^. The pfiysical radius of tlie fuzzy 
is 



To 



(72) 



Similar to tlie D3_LD1 system in which the wormhole-hke solutions connecting a D3 and an 
anti-D3 exist, there would be parallel M5-anti M5 solutions in the ABJM model. The equation 
nl) can be rewritten as 



Cl, 



(73) 



where Ci is an integration constant. Let us consider Ci = — a/g, /o > case. In this case, 
because (/')^ = ap — af^ > 0, the solution is valid only in the region / > /q. Considering 
minus sign in (1731) . / is a monotonically decreasing function of X2- Remarkably, /' = at 
/ = /o that happens within a finite value of X2. The equation (1731) can be integrated to give 



X2 



+ 



df 



f 



In terms of the physical radius, X2 can be rewritten as 



X2 



+ 



2(A^ 



^/aT2 



dR 



JTLq 



(Rl 



2{N 



To 



(74) 



(75) 



At large i?, this solution is approximated by the BPS fuzzy funnel solution, namely, an M5- 
brane at X2 = Xoo- At R = Rq > 0, {x2 = X2{Ro)), the fuzzy funnel ceases to decreasing 
its radius. As in the case of Dp/anti-Dp [18], this solution can continue past this point by 
introducing another part of the solution 



X2 = Xoo + 2AX2 



2(A^ 



y/aT2 



dR 



R 



^/R^ 



Rl 



AX2 = X2{Ro 



(76) 



The solution ( 1761) is smoothly connected to the solution ( 1751) at R = Rq. The solution ( 1761) is 
the solution of (173!) for plus sign so that R increases from Rq. For large R, the solution (!76|) 
is approximated by another 5-brane at X2 = x^o + 2Ax2. As in the case of the D3_LD1, this 
another 5-brane can be interpreted as an anti M5-brane. This is because the orientation of 
the brane smoothly becomes opposite. See fig.[7]for a schematic picture. The integral in the 
solution is evaluated as 



dR 



Rl 



QRl 



6Rl 



dy y 6 (y 



dy y ^{y- l)' 



2' 3 



B 



R^-R^n 1 1 



i?6 



'2'3 



(77) 



where we have changed the integration variable from R to y = Rq ^R^. Be careful that there is 
no point X2 corresponding to i? = 0. From this expression, we have 

. k 1 



AX2 = X2{Ro) 



Xr 



2{N 



27rT2 6Rl 



B 



1 1 

2' 3 



OC Rn 
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Figure 7: A schematic picture of the solution (^^, (17U|) . 



Note that this solution is not BPS but for Rq = limit, one of the two 5-branes goes to infinity 
and we recover the BPS fuzzy funnel solution which is nothing but an (anti) M5-brane. For 
large separation of M5/anti-M5, Ax2 — > oo, the throat shrinks to zero-size. This is similar to 
fundamental strings between a D-brane and an anti-D-brane. 

This fact can be confirmed from the M5-brane point of view. A scalar field X in a single 
jo-brane world-volume is governed by the Nambu-Goto action 



Sp = Tj <F+^i + d.Xd^^X. (79) 



Here X describes fluctuation along a transverse direction to the p-brane and Tp is a tension of 
the p-brane. The static, spherically symmetric equation of motion of X is 

= 0, (80) 




where r = y + ■ ■ ■ is a radial coordinate in the p-brane world- volume and the prime in the 

equation stands for the differentiation with respect to r. The equation ( IHOl) can be integrated 
to give 



p- 

X{r) = I dr- 



,2p-2 ^2p-2 
'O 



Here is an integration constant. This precisely match with the solution (!75|) with p = A. Note 
that because the M5-brane shares the one-dimension with the M2-branes, p = 5 is effectively 
reduced to p = 4. 

The energy of the solution is evaluated as 

E = / dtdx^ — dR , ° =. (82) 

47rAr-i7 \k ^R^ - Rl ^ ^ 
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Figure 8: An analytic plot of X2(/) with k = fo = 1. The solid and dashed lines correspond to 
plus and minus signs respectively in the equation (1731) . 

For large and large R region, the energy is approximately given by 



brane. Note that the factor k in the volume appears due to the orbifolding Z^. This is another 
evidence that the 5-branes located at the divergent points can be interpreted as an (anti) M5- 
brane. To clarify whether the brane we are considering is M5 or anti-M5, we need to evaluate its 
charge associated with 3-form in the eleven- dimensional supergravity. Indeed, for the D-string 
case, we can show that the configuration of the double funnel becomes the source of the R-R 
4-form potential and the opposite world-volume directions correspond to different signs of the 
charge [21]. Unfortunately, in our case, we do not know the general coupling of 3-form in the 
multiple M2-brane world-volume theory and in order to clarify the issue more explicitly we 
need further study of this coupling to evaluate the charge corresponding to this configuration. 
Currently, we are unable to determine this charge and this issue is beyond the scope of this 
paper. 

4.2 Cusp solution 

Let us consider ci = a/g, fo > case in the equation flTSj) . In this case, the fuzzy funnel 
collapses at finite X2 (the minus sign solution in ( 173|) ) where the fuzzy funnel is pinched off. An 
analytic profile of the solution is plotted in fig. [HI 

As discussed in [19,21] for D3_LD1 case, there is no singularity at i? = and the solution can 
pass through this point smoothly provided the negative radius has physical meaning. Alterna- 
tively, once we want to keep the radius positive, the solution can be extended by introducing 
the solution for plus sign in ( 1731) . This is a cusp solution which would represent M2-branes 
connecting an M5-brane and another M5-brane. This is because the orientation of the brane 
world- volume does not change. A typical profile of this solution is illustrated in fig.O This 
situation is analogous to the one discussed in [19] in the D3-brane context. Especially, for 



E = T5Vol(M5), 



(83) 



where = and vol(M5) = J dx 




R^dR are the tension and the volume of an M5- 
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Figure 9: A schematic picture of the cusp solution. 



R> region, the integral can be evaluated as 
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R ^fWTM 
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Contrast to the M5/anti-M5 solution, there is a point X2 corresponding to = 0, - a shrinking 
point of fuzzy S^. The distance between two parallel M5-branes is given by 



Ax2 = 2{N -1] 



k 



27rT2 6Rl 



B 



1 1 

6' 3 
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1 1 1 

2' 6' 3 



OC Rn 



^5) 



A dual picture from the viewpoint of the M5-brane would be possible. As in the case of the 
D3-brane analysis, turning on the self-dual gauge field in the M5-brane world-volume [32, 33] 
would play an important role to construct the cusp solution presented here. Although there 
are two parallel M5-branes, this solution is not BPS. The BPS configuration is recovered in the 
limit Ro ^ where one of the M5-brane goes to infinity. 



4.3 Non-BPS solutions in mass deformed ABJM model 

Let us briefly discuss the static solution of the equation of motion in the mass deformed ABJM 
model. With the ansatz = f{x2)S^, / G R, the equation of motion (1161) can be rewritten 
as 

There is an integration constant c defined by 



c=^f"-V{f). (87) 
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Note that for c = 0, this equation reduces to the BPS equation 



f = ±ff{f-v'), (88) 

and the solutions are given by f± in ffTHj) . As in the case of the time-dependent solution, it is 
convenient to rewrite the equation of motion with respect to t = p 

= 8f[Vif) + c] = [r(r - + cv'] , (89) 
This equation can be simplified by using the following ansatz 

t(x) = — — r- . 90 



With this ansatz, the equation fl89|) reduces to the equation for the function h 

ih'f = Ah'-g2h-gs, (91) 
where the parameters g2 and gs are given by 

^2 = ^(6c+l), g,^-A(^c' + ^c+^y (92) 

From equation 0911) . we find that the function h is given by the Weierstrass's elliptic function 
p. Therefore, the static solution of the equation of motion is 

r = f = . (93) 

3p{m{x — Xq)) — 1 

For the values of the parameter c < and c > 0, the solution has similar profiles as fig. [7] and 
fig. [H] respectively. By taking the limit c — with adjusting the parameter xq appropriately, 
we can reproduce the BPS solutions f± in (|TH|) . 



5 Conclusion and discussions 

In this paper, we have investigated several classical solutions of the AB JM model. The equation 
of motion for the scalar fields can be solved by the help of the BPS matrices. Due to the 
noncommutative nature of the matrices, the solutions generically exhibit fuzzy structures. 

In the first part of the main text, we studied time evolutions of fuzzy spheres particularly 
focusing on fuzzy S^s. We found analytic solutions both for the massless and massive cases. 
For the massless case, there is an oscillating fuzzy sphere solution which collapses into zero size 
then expands again. The decay time is analytically derived. For the massive case, there are 
solutions having different behaviors according to the value of the mass and the initial condition. 
Remarkably, if the mass has an appropriate value, there are solutions that do not decay into 
zero size but keep its radius being positive. This is different from the massless case where such 
solutions do not exist. 
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In the second half of the main body, we investigated purely spatial, static, generically non- 
BPS solutions. Similar to the D3_LD1 system, we found two kinds of solutions. One is the 
wormhole-like solution with its cross-section fuzzy S^. This solution is given by smoothly 
connecting two fuzzy funnels. The two fuzzy funnels expand into two M5-branes with opposite 
world-volume direction. Therefore one of the M5-branes can be interpreted as an anti M5- 
brane. The solution can be also analyzed from the viewpoint of dual M5-brane world-volume 
picture. The other is the cusp solution connecting two 5-branes. In this solution, the fuzzy 

is pinched off at a point and then growing into the other 5-brane. Because the 5-brane 
world-volume direction does not change, the two 5-branes are interpreted as two M5-branes. 

There are a few comments on the study performed in this paper. The solutions examined 
in this paper are solutions of the effective action at leading order in the derivative expansion 
of a non-linear action. This corresponds to the leading order Yang-Mills part in the expansion 
of the DBI action. Actually, based on the novel higgs mechanism first proposed in [34], when 
a scalar field develops a VEV and taking the large VEV and large k limit with fixed v/k, the 
ABJM model reduces to the three-dimensional jV = 8 U{N) super Yang-Mills action - the 
leading order of the effective action of N D2-branes in type IIA string theory [13,35]. In the 
discussion of D-brane world-volume effective action, such a leading order solution (including 
the wormhole and cusp solutions) frequently gives a solution in the full non-linear level and 
the solutions found in this paper would keep being solutions at the full non-linear M2-brane 
effective action. Some non-linear version of multiple M2-brane effective theories were proposed 
in [36-38]. 

Another important point is to find M2-brane effective action in general supergravity back- 
ground. The ABJM model is interpreted as an multiple M2-brane effective action in flat 
space-time C^/Zfc, and couplings with the background supergravity field cannot be described 
correctly in this model. It is important to find the correct coupling of these supergravity fields 
in the multiple M2-branes. 

In this paper, we analyzed the simplest class of the solutions, ^ 0, Wa = = = 0. 
It would be interesting to investigate the effects of non-zero gauge fields like the situation 
studied in [21] where (p, g)-strings attached to a D3-brane. 
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A BPS matrix 



The BPS matrices 



(A = 1,2) satisfy the following relations, 

oA oB o\ oA oA o\ oB 




(94) 
(95) 
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The explicit form of the matrices was first found in [26] . This is given by 

{Sl)mn = Vm- 16m.n, ('S'|)mn = VN - m6m+l,n, (96) 

S^Sl = diag(0, 1, 2, ■ ■ ■ , iV - 1) = SlS\ (97) 

S^Sl = diag(0, iV - 1, iV - 2, ■ ■ ■ , 1), (98) 

5|52 = diag(A^-l,A^-2,--- ,1,0), (99) 

S'^S\ = diag{0,N,--- ,N). (100) 

From this expression, we have 

Tr^^^;!^ = Tr^;^^^ = N{N - 1). (101) 

Let us see the hnearization of the matrices S"^ in the right hand side of the massless equation 
of motion (fT2l) . Assuming an ansatz = f{x)S'^, the matrix structure on the right hand side 
of ([121) is 



3(5^4)^5^ + 35^(45^)' - 25^4(5^4)5^ - 2(5^5^)5^(51^5^) - 25^(5^5^)51^5^ = M^. 

(102) 

The BPS matrices satisfy the relation 

{S^'SlfS'^ = {S''Sl)S'^ + {S''Sl,)S'^{SlS^), (103) 
S'^iSlS'')^ = -S'^{Sl,S'') + {S''Sl)S'^{SlS'='). (104) 

Therefore 

3(5^51^)^5^ + 35^(5lj5^)2 = 35^ + 6(5^51^)5^(5^,5^). (105) 

Then 

= 35^ + 4(5^5lj)5^(5j5^) - 2S'^Sl{S^ S^^^)S'' - 25^(5^,5^)51^5^. (106) 
Using the relation 

5lj5^ = SUS^Sl)S'' -{Slsy, (107) 
5^5lj = (5^5j)2 - 5^(5j,5^)5lj, (108) 

we have 

4(5^5^)5^(5^,5^) - 25^5^ (5^5^"^^) 5^ - 25^(5j,5 C)5lj5^ 
= 4(5^5lj)5^(5j,5^) - 25^(5lj5^) + 2(5^51^)5^ + 25^(51^5^) 
-2(5^5lj)5^(5j;5^) - 2(5^5lj)5^ - 2{S'' Sl,)S'^{SlS^) 
= 0. (109) 



Finally, we find 



M^ = 35^. (110) 



For the massive case, there is only an additional term coming from the second term in the right 
hand side in fllSj) . This can be trivially linearized by using the relation 
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B Moving fuzzy funnels and a domain wall 

As discussed in section 1, the equation of motion is world- volume Lorentz invariant. Therefore 
one way to obtain a non-trivial time-dependent solution is Lorentz boosting of known purely 
spatial solutions. In this appendix, we show that moving fuzzy funnels and domain wall so- 
lutions with constant velocity are obtained by boosting known BPS solutions. Consider an 
ansatz 

= /GR, (111) 
then the equation of motion for massless case ( fT2i) reduces to 

The BPS solution 

/(x)=/(x2) = yj(x2 + Xo)-^ (113) 

satisfies the equation (fT2ll . In addition to the static solution, we can explicitly show the following 
solution satisfies the equation, 



f{x2,t) = ^—(jX2±Vl'-U) \ (114) 

Here < 7 < 1 is a real parameter. This is just the Lorentz boost of the BPS solution flllSp 
with 

x''^ = -f(x^ -vt), 7= , ^ =. (115) 



^ /^2 _ I 

This solution would represent moving an M5-brane with constant velocity v = ±c-^^^ — . c is 
the speed of light. Let us check this fact. The cross-section of the solution f lll4p at fixed X2,t 
is a fuzzy with radius 

fl^ = (116) 

27rT2 ax2 ± va^ — It 

Then the energy of the configuration (11141) is evaluated as 



E = I dtdx^dx^Ti 



JL^ [ dtdx^ ( R^dR. (117) 



27iN ~1 J \ k 

At large- A^, this reduces to 



E = ^=== {N^oo), (118) 
a/1 — v^/c^ 
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where 

M = T5Vol(M5) (119) 

is the mass of an M5-brane. This energy precisely match with the one of an M5-brane moving 
with velocity v < c. Similar to the massless case, we obtain moving deformed fuzzy funnel and 
domain wall by Lorentz boosting known massive BPS solutions 

h{x„t) = ^{l±e-'-'^y\ (120) 
where x' is given by the equation flllSp . 



C Time-dependent solution in a simple form 

In this section, we derive a simple form of the time-dependent solution which is valid for 
arbitrary values of /q. The equation of motion (l32l) can be rewritten as 



This is nothing but the equation of motion for the particle in the potential V{f). Therefore we 
can define the conserved energy as 

E^^f' + V{f). (122) 

If we assume that /(to) = /o and /(to) = 0, the energy is given by 

E = V{fo). (123) 

The equation (f^ = —8p[V{f) — V{fo)]) can be cast into a simple form by using the 
following ansatz 



(/o^-^^)(3/o^-^^) ,_1 , 2/oW 



V 



V b + h{m{t - to)) J '3 ^ ' 

Then, the equation for the function h is given by 

= Ah' - g2h - g,, (125) 

where the parameters g2 and g^ are given by 

g2^-li2c-l), g,^±{3c^-Qc + 2), c^^JM_^, (i26) 

The solution to the equation fll25p is nothing but the equation for the Weierstrass's elliptic 
function p. Therefore, the time-dependent solution for r is given by 



- - ^» I 1 (127) 



This solution is valid for arbitrary value of /q. 
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